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Abstract. We study the relaxation dynamics of the one-dimensional Tomonaga- 
Luttinger model after an interaction quench paying particular attention to the 
momentum dependence of the two-particle interaction. Several potentials of different 
analytical form are investigated all leading to universal Luttinger liquid physics in 
C^ , equilibrium. The steady-state fermionic momentum distribution shows universal 

behavior in the sense of the Luttinger liquid phenomenology. For generic regular 

potentials the large time decay of the momentum distribution function towards the 

^ ] steady-state value is characterized by a power law with a universal exponent which 

O ' only depends on the potential at zero momentum transfer. A commonly employed 

ad hoc procedure fails to give this exponent. Besides quenches from zero to positive 
interactions we also consider abrupt changes of the interaction between two arbitrary 
K^ values. Additionally, we discuss the appearance of a factor of two between the steady- 

p^ , state momentum distribution function and the one obtained in equilibrium at equal 

CO ' two-particle interaction. 
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1. Introduction 

The experimental progress in controlling and manipulating cold atomic gases [1], which 
in certain parameter regimes form strongly correlated quantum many-body systems, led 
to tremendous activities with the final goal of theoretically understanding the quench 
dynamics in such systems (for a recent review, see Ref. [2]). In a sudden quench at 
least one parameter of a given Hamiltonian is switched abruptly at time t = — the 
Hamiltonian before the quench is Hi, the one afterwards H{. At t = the system is 
assumed to be in the ground state (we consider temperature T = 0) of Hi and the time 
evolution for t > is performed with respect to H{. We here study 'global' quenches of 
the two-particle interaction. 

Describing the time evolution of a correlated quantum system out of a 
nonequilibrium state poses a formidable challenge. It is reasonable to consider the case of 
one-dimensional (Id) systems first, as in Id a variety of analytical as well as numerical 
methods exist which allow for controlled access to equilibrium correlation effects in 
specific models [3l HI |5l Ej . Several of those techniques were recently extended to study 
the nonequilibrium problem at hand [3, [HI [HI [10], [TT], [12l [131 [HI [ISl [IE]. Furthermore, 
in Id chains virtually all many-body states of matter of current interest, such as 
Mott insulators, superfiuids, superconductors and charge- as well as spin-density-wave 
states can be realized. Quenching between these states is expected to be particularly 
interesting. 

We here focus on fermionic systems. Even if the system stays metallic in the 
presence of two-particle interactions the Fermi liquid concept breaks down. Instead 
a wide class of Id models shows Luttinger liquid (LL) behavior on low-energy scales 
[3],[11[5]. Equilibrium LL physics is characterized by the 'universal' power-law decay of 
certain correlation functions in space-time with exponents which for spinless and spin- 
rotational invariant models can all be expressed in terms of a single number K. This 
LL parameter in turn is a function of the microscopic details of the model considered, 
in particular the strength of the two-particle interaction. For noninteracting fermions 
K = 1 and < i^ < 1 for repulsive ones; the case we focus on. 

The Tomonaga-Luttinger (TL) model is the exactly solvable effective low-energy 
fixed point model of the LL universality class [T71 [T8l [31 [U [5]. It thus plays a similar 
role as the free Fermi gas in Fermi liquid theory. The model has two strictly linear 
branches of right- and left-moving fermions and two-particle scattering is restricted to 
processes with small momentum transfer |g| <^ kp, with the Fermi momentum kp. These 
processes as well as the kinetic energy can be written as quadratic forms of the densities 
of right- and left-moving fermions which obey bosonic commutation relations. In most 
calculations the momentum dependence of the low-momentum scattering processes g2 
and 5^4 (in the so-called g-ology classification [18]; see below) are (partially) neglected 
and momentum integrals are regularized in the ultraviolet introducing a convenient 
cutoff 'by hand'. 

In the bosonization procedure the interacting Fermi system is mapped onto a model 
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of free bosons. The time evolution of the latter is trivial and the dynamics of correlation 
functions of the fermionic densities (at small momenta) after an interaction quench 
can be accessed directly. In addition, the fermionic field operator can be written as a 
(highly nonlinear) function of the bosonic eigenmodes an [HI [21 HI [5] such that the 
time evolution of fermionic correlation functions can be computed exactly as well. The 
TL model thus constitutes an ideal play ground for studying the dynamics resulting 
from an interaction quench. 

The time evolution of the fermionic single-particle Green function Gt{x) = 
(^ip^{x)ip{0)'^ , > as well as the density correlation function (at small momenta) after 
suddenly turning on the interaction in the spinless TL model was first studied in Ref. |19j . 
Here ^^'^''(x) denotes the field operator and (. . .) ,^^ the expectation value with respect 
to the time-dependent density matrix p(t). In Ref. [19] it was found that at large t 
the Green function approaches a time-independent stationary limit. At T = the 
stationary Green function shows power-law behavior as a function of the position x 
with a K dependent exponent. Power-law decay of \im.t^ooGt{x) = G^^{x) at large \x\ 
translates into a typical LL power-law behavior of the stationary fermionic momentum 
distribution function close to kp, \n^^{k) — 1/2| r^ \k — /cfT'S with 7st being a function of 
K and thus depending on the strength of the interaction. Interestingly, the steady-state 
exponent 7st differs from the one of the ground-state momentum distribution function 
7gs at the same interaction strength. For finite times n{k, t) has a Fermi liquid-like 
jump at kp, with a Z-factor which vanishes as a power law in t: Z ~ t"'''"'. Further 
aspects of the quench dynamics of the TL model or closely related ones were discussed 
in Refs. [201 [2ll [221 [231 [21] . 

In Ref. [19] an ad hoc ultraviolet regular izat ion was used. It is widely believed that 
(partially) neglecting the momentum dependence of the interaction and regularizing 
momentum integrals as convenient has no effect on the low-energy equilibrium physics 
of the TL model. This is indeed correct if all energy scales are sent to zero |25j . 
The such obtained results for the dependence of power-law exponents of specific 
correlation functions on the LL parameter K become universal and are valid for 
all models falling into the LL universality class. It is however questionable if the 
same reasoning holds when considering quenches. High energy processes and thus 
the full momentum dependence of the two-particle interaction might matter in this 
nonequilibrium situation. To investigate this issue we keep the momentum dependence 
of the g2 and g4 processes — rendering any ad hoc ultraviolet regularization superfluous 
[25] — and consider interaction potentials of different analytical form. We exactly 
compute Gt{x) and its Fourier transform, the fermionic momentum distribution function 
n{k, t), of the spinless TL model after a quench out of the noninteracting ground state. 
We first show that independent of the details of the momentum dependence of the 
interaction the steady-state momentum distribution as a function of k is characterized 
by a power-law nonanalyticity at kp with the exponent '-fst{K). The LL parameter K is 
a function of the potential at zero momentum transfer only. Therefore n^^{k) of the TL 
model is universal in the LL sense. Similarly, the Z-factor shows a universal power-law 
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decay in time. We then proceed and show that for generic regular potentials also the 
asymptotic time dependence of n{k, t) has universal aspects. For fixed k ^ kp we find a 
power-law decay towards the steady-state value as a function of t with the /c-independent 
exponent l-|-7st- The ad hoc procedure instead gives the exponent l-|-37st/2 and is thus 
insufficient for studies of the generic dynamics. The power-law decay is overlaid by an 
oscillation with a frequency which depends on the momentum k considered as well as 
the momentum dependence of the potential. In addition, we investigate a box shaped 
(in momentum space) potential and show that the asymptotic dynamics is dominated 
by its discontinuity leading to a power-law decay with exponent 1 (independent of the 
interaction strenght). We raise the question if the universality in the steady state as 
well as the large time dynamics of the TL model after an interaction quench extend to 
other models of the (equilibrium) LL universality class. Besides quenches from zero to 
positive interactions we briefly consider quenches between two interactions of arbitrary 
strength. 

Our paper is organized as follows. In Sec. [2]we present the model and give details on 
how to compute the single-particle Green function after a quench. Particular emphasis is 
put on the ultraviolet regularization by momentum dependent two-particle potentials. 
Different model potentials are introduced. In Sec. |3] we compute and compare the 
steady-state and equilibrium Green and momentum distribution functions. For small 
interactions they differ by a factor of two which was earlier discussed in the context of 
pre-thermalized states [261 1271 [28]. In Sec. Hlwe present our results for n{k, t) and show 
that they partly depend on the form of the potential considered. We also discuss results 
for the standard ad hoc regularization. Finally, our findings are summarized in Sec. [5l 

2. Model and methods 

2.1. The Tomonaga-Luttinger model and bosonization 

The interaction energy of a translational invariant system (periodic boundary 
conditions) of Id spinless fermions interacting via a potential which only depends on 
the distance between the two scattering particles can be written as 

1 Oyr 

-^Pot = 7 X^ v{qn)PnP-n + hi{N) , Qn = — U , U E Z , (l) 

n.>0 

with the density operator 

m 

and the Fourier transform v{k) of the two-particle potential V{x). Here the Cn denote 
fermionic momentum space annihilation (creation) operators, hi{N) contains terms 
which depend on the particle number operator A^, and L is the length of the system. 
After linearization of the single-particle dispersion around the two Fermi points at ±/cf 
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the kinetic energy reads 

-f^kin = "^Vpikn - /CF)ci,+C„^+ + ^{-VF){kn + ^f)c1,_C„_ + h2{N) , (3) 

n>0 n<0 

where we already introduced independent right- (c„>o,+) and left-moving (c„<o,_) 
fermions [5]. The Fermi velocity is denoted by vp. In the next step one supplements the 
Hilbert space of the right movers by states with negative momenta and the one of left 
movers by states with positive momenta. The linearization and addition of states does 
not change the equilibrium low-energy physics. From now on we drop terms containing 
A^ as they are irrelevant for our considerations. The kinetic energy H]^in can then be 
written as a quadratic form in the densities pn,± of the right and left movers defined in 
analogy to Eq. ([2]) [5]. With 



Pn, 
Pn, 



for n > 
for n < 



(4) 



one obtains 



The bn obey the standard bosonic commutation relations. Replacing p„ — )■ /?„ + -|- p. 
in Eq. ([T]), using Eq. (^ we obtain the Hamiltonian of the TL model 

i^TL = E W U + ^) {bib^ + bij_:) + k:^ (bibi 

n>0 L V / 



(5) 



b-nbn 



(6) 



Distinguishing between intra- and inter-branch scattering processes one often replaces 
the potential v{k) in the first term by a function g^i^k) and the one in the second by 
an independent function g2{k) [H]. For simplicity we here refrain from doing so and 
assume g2{k) = g^ik) = v{k). This has no effect on our main results. 



2.2. Eigenniodes and dynamics 

By a Bogoliubov transformation the 'bosonized' Hamiltonian i^TL can straightforwardly 
be diagonalized 

H^h = Y^ uj{kn) aian + ^gs (7) 

by introducing the eigenmodes 

a„ = c{kn)b^, - s{k„)bl„ ^ h^ = c{k„)a^ + s(/c„)aL„ (8) 

with 



S [kn 



l+t)(fc„)/2 



- 1 



c2(fc„) = l + s2(A:0 



L>j{k^ 



i) = Vf \kn\ >yi + v{kn) , V{kn) = — ■ 



TTVf 



(9) 
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In the noninteracting limit s'^{kn) -^ and uj{kn) -^ fpl^nl- For physical reasons the 
Fourier transform v{q) of the two-particle potential must vanish on a characteristic scale 
denoted by gc- This implies 

lim s^{k) = , lim c'^{k) = 1 , lim ^^ = vp . (10) 

fc— >00 k^fOD fc— >oo \k\ 

The LL parameter K of the TL model is 

K = [i+v{o)r'^^ . (11) 

It thus only depends on the (dimensionless; see Eq. ([9])) potential at momentum g = 0. 
This shows that the TL model is a LL with < A' < 1 only if v{0) > 0; we here focus 
on two-particle potentials with this property. Using Eq. iQ one finds 

,2(0) = 1 (K + ;^-i _ 2) , c\0) = 1{K + K-' + 2) . (12) 

It turns out to be useful to introduce a renormalized velocity vp and its dimensionless 
analog vp 



dujik) 
v-p ^ 



, ^F = -. (13) 

fc=0 ^F 



dk 
Using Eq. (I9l) we find 

vp = Vl + C'(0) . (14) 

The ground state of i^TL is given by the vacuum with respect to the eigenmodes 

I vac (a)) and 

Egs = -2t;F Y^ kn S^{kn) ^/iTWn) (15) 

n>0 

is the ground-state energy. For vanishing two-particle interaction E^^ = 0. 

The time evolution with respect to H^ = H-pi^ of the eigenmode annihilation and 
creation operators in the Heisenberg picture is now trivially given by 

a„(t) = e-'"('=")*a„, ai{t) = e'^^'^-^'al (16) 

and the one of the hn can directly be obtained using this and Eq. (|8]). 

2.3. Bosonization of the field operator and time evolution of the Green function 

To obtain expectation values of fermionic operators, such as the momentum distribution 
function of the right movers (the left movers can be treated similarly) 

n{kn,t)= dxe^'-^Gtix), Gi(x) = (v^i(x)^+(0) ) , (17) 

J-L/2 \ ' P(*) 

one has to 'bosonize' the fermionic field operator of the right-moving particles 



*'iw = 4fE 



X/CtiX T 



Vl^^ c,,+ • 
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^—ixir/L 



^l(x) = ^-i^H-)uie-"^i-) ^ (19) 



with 



where [/^ denotes a unitary fermionic raising operator which commutes with the hil and 
maps the A^-electron ground state to the (A^ + l)-electron one. In the computation of the 
Green function the fermionic operators lead to the phase factor exp {—ik-px) appearing 
below (see Refs. [5] and [25] for details). 

With the initial density matrix p\ = |vac(6)) (vac(&)| corresponding to the 
noninteracting ground state and the time evolution given by the interacting Hamiltonian 
Eq. ([7]) we obtain 

Gt{x) = (iji{x)i^+{0)) = (vac(6)|^^(x,t)^H_(0,t) |vac(6)) 

E ^"'^^"^"'^^"^ (cos [Kx] - 1) (1 - cos [2u{K)t])\ , (21) 
l~e^^-) L^ ^ J 

where we used Eqs. (IHl) and (IT6l) as well as the Baker-Hausdorff relation. To prevent 
recurrence effects at large times we take the thermodynamic limit 

Gt{x) = — ^^exp ( Hdk ^^ ^^^"^ ^^^ (cos [kx] - 1) (1 - cos [2uj{k)t])\ . (22) 

We stress that because of Eq. (TTOl) the momentum integral is convergent at large k and 
does not require any regularization. 

A more general situation arises if one starts at t = in the ground state with two- 
particle potential Vi{q) and performs the time evolution in the presence of the potential 
V{{q) (an analogous situation for bosonic LLs is discussed in Ref. [29]). Applying the two 
Bogoliubov transformations from the noninteracting eigenmodes to the ones with Vi{q) 
and ff(g) given by relations analogous to Eqs. (IHl) and ([9]) one can generalize Eq. (l22l) 
to (in self-explaining notation) 

t e-^^- r r cos[kx]-l 

Gt [x) = — exp <.2dk 




2% X + iO I Jo k 

X I ciik) [cf{k)s,{k) - Si{k)c,{k)f + s^.ik) [ci{k)c,{k) - s,{k)s,{k)y 



+2sf{k)cf{k) [c{{k)si{k) - S{{k)ci{k)] [cf{k)ci{k) - S{{k)si{k)] cos [2uj{{k)t] J \ . (23) 

For V[{q) = 0, implying Ci{k) = 1 and Si{k) = 0, this expression reduces to Eq. ( l22l) . 

Equation (l23l) also covers the interesting situation in which one considers the ground 
state with nonvanishing Vi{q) as the initial state and performs the time evolution with 
the noninteracting Hamiltonian, that is for V{{q) = and thus Ci{k) = 1 and Sf{k) = 
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In this case the fermionic momentuin occupancy is conserved. Accordingly the 
time dependence in Eq. f l23|) drops out and the Green function reads 

Gt{x) = f ^exp ( rdk2sKk) S^^M^I = G^x) . (24) 

It corresponds to the equilibrium ground-state Green function G^^{x) of the TL model 
at interaction Vi{q). Although the time evolution is with a noninteracting Hamiltonian 
the system is characterized by the LL Green (and momentum distribution) function of 
the initial state for all times. 

2.4- Momentum integrals and potentials 

Computing expectation values in the TL model in and out of equilibrium one 
regularly faces momentum integrals of the type appearing in Eqs. (l22l) . ( l23l) . and (12^ 
[31 m El [251 [ini [H]- To evaluate those one often assumes that the prefactor of the 
trigonometric functions has a convenient exponential form. In the case of the time 
evolution (under the interacting Hamiltonian) out of the noninteracting ground state 
one sets (see Eq. (|22Dll|| 

As\q)c\q) = (^^e"!^!/"^ . (25) 

Using Eq. (Q one can solve for the momentum dependence of the underlying two-particle 
potential and obtains 

i,(q) = 2^2e-lg|/gc ^ 2 x/gV^HTgT yf^^V^kCV . (26) 

The momentum-space potential decays exponentially for \q\/qc ^ 1- Even with this 
special choice of the momentum dependence of the two-particle potential the integral in 
Eq. ( I22l) cannot be solved analytically as uj{k) appears in the argument of the cosine. 
One therefore linearizes the dispersion and makes the replacement 

uj{k) -^VF\k\ VA;, (27) 

with vp given in Eq. (fT3l) . Equations (!26l) and ([27]) form one of the possible ad hoc 
procedures used in the literature [31 [H [5l [251 [IH] to make analytical progress; others are 
applied as well. With these a closed analytical expression for Eq. ( l22l) can be given (see 
Sec. [1]). It is usually believed that this replacement (and similar ones) does not change 
the equilibrium low-energy physics, which is indeed correct if all energy scales are sent 
to zero [25]. However, in nonequilibrium the dynamics is affected by the high-energy 
modes and thus the full momentum dependence of the potential becomes important. 
The replacement was still used in Refs. [T9l [2^ based on the expectation that at least 
in the weak coupling limit it will not significantly affect the time dependence. Here we 
do not rely on this approximation and exactly evaluate the integral in Eq. (l22l) for the 
potentials 

/ \ [ 1^ for \q\ < qc ,^^. 

^hoM = S n f k ' 28 

tor \q\ > qc 



I In ground-state calculations one instead sets 2sf{q) = g^ exp (— Igl/^c) [30], see Eq. ([21 
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Vuss(g)=^^e-('^/^=)V^ (29) 

^^exp(g) =t;e-l''l/'^% (30) 

V 
^^quart(g) = "—4 ■ (31) 

l + (t) 

All potentials have the same g = value v{0) = v and thus the same LL parameter K 
and renormalized Fermi velocity -Op (see Eqs. ( ITTl) and dHD). In equilibrium they give 
the same low-energy LL physics. To obtain this also for the potential Eq. (126|1 one has 

to choose 

1 V , ^ 

In Sec. m we show that the quench dynamics has aspects which are equal for different 
choices of the potential while other features depend on the potential. Before discussing 
this, in the next section we compute G^^{x). The steady-state Green function and 
thus n^^{k) turns out to be universal in the sense of the LL phenomenology. We 
furthermore compare n^^{k) to the corresponding equilibrium momentum distribution 
function computed for the same interaction strength. 

3. Steady-state and equilibrium expectation values 

3.1. Infinite-time lim,it and equilibrium 



The infinite-time steady-state value of Gt{x) Eqs. (1221) and ( 1231) can be obtained 
straightforwardly. We first consider the quench out of the noninteracting ground state. 
For fixed large t the cosine term in Eq. ( 122|) with argument linear in t becomes a rapidly 
oscillating function of k and averages out. In the limit t — )■ oo we thus end up with 

G'\x) = lim Gt{x) = — ^^ exp | Hdk ^^ ^^^"^ ^^^ (cos [kx] -1)1. (33) 



t^oo 2-11 X + iO iJq k 

The fermionic momentum distribution function in the steady state n'^^lk), that 
is the Fourier transform of Eq. fl33l) . can close to kp be computed without any 
specific assumptions for the two-particle potential using asymptotic analysis [31]. One 
can closely follow the steps of Ref. [25] for the equilibrium ground-state momentum 
distribution function obtained from Fourier transforming Eq. (124|) . Next we briefly 
outline those. The leading large |x| behavior of Eq. (^^ is given by the integrand at 
small momenta and dominates the Fourier transform close to kp leading to 

\n^'{k) - 1/2| ~ AA; (34) 

independent of the details of the potential [25]. Here we have introduced the relative 
dimensionless momentum 

Ak = '^^^ . (35) 

The equilibrium anomalous dimension 7gs reads (we drop the index i) 

^^, = 2s'{0) = l{K + K-'-2) (36) 
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Equation fl34p holds for all two-particle potentials within the TL model as long as 
< 7gs < 1; that is for sufficiently small interactions. It implies a power-law singularity 
of the first derivative of n^^{k) at kp. For larger interactions n^^{k) goes linearly through 
kp and singularities appear in higher order derivatives [25] . 

In complete analogy one obtains for the steady-state momentum distribution 
function rf^(k) after a quench out of the noninteracting state 

\n'\k) - 1/2| ~ Ak ^" (37) 

with the nonequilibrium anomalous dimension 



7st = 4.^(0)c^(0) = - 



4 1 + v{0) 
= - [K^ + A'"2 _ 2) . (38) 

For Eq. ( 137|) describing the leading behavior close to kp the same restrictions as outlined 
in connection with Eq. ( 134|) hold (with 7gs replaced by 7st). In the specific case of the 
potential Eq. ( I26l) the results Eqs. (1371) and (l38l) were obtained in Ref. [19]. The ad 
hoc replacement Eq. fl27|) is not required here since in the limit t — )■ oo, u{k) drops out. 
Interestingly, 7st in the steady state of the TL model differs from the exponent 7gs found 
in the ground state of -ffxL at the same interaction strength. In the next subsection we 
return to this issue. 

Quenching between two repulsive interactions of arbitrary strength the steady-state 
Green function follows from Eq. fl23|) by dropping the cosine term with time dependent 
argument 



271 X + iO I Jo ^ 

X (cUk) [c,{k)s,{k) - Si{kMk)f + sl{k) [ciikMk) - s,{k)s,{k)f) \ . (39) 



Independently of the details of the momentum dependence of the initial and final 
potential Fourier transforming again leads to Eq. ( 137|) with the nonequilibrium 
anomalous dimension given by 

7st = 2 (c2(0) [cf(0)si(0) - Sf(0)ci(0)]' + 5^(0) [cf(0)ci(0) - Sf(0)si(0)]') . (40) 

Using Eq. flT^ (supplemented by the indices i and f ) 7st can be written as a function of 
the LL parameters Ki and K[ associated to the two interactions. As this does not have 
a simple form we refrain from presenting it here. 

Power-law behavior of correlation functions is a typical feature of LL physics. 
Within the TL model we thus find universal behavior known from equilibrium: the 
steady-state momentum distribution function close to kp is independently of the details 
of the momentum dependence of the potential characterized by a power law with an 
exponent which can be expressed in terms of K or K\ and Ki. It would be very 
interesting to investigate if this universality of the steady-state expectation value extends 
beyond the TL model. For this it would be necessary to analytically or numerically 
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compute n^^{k) for other models from the LL universahty class, e.g. lattice models 
such as the model of spinless fermions with nearest-neighbor hopping and interaction, 
for which K is known from other considerations (e.g. the Bethe ansatz or numerics) 
O m [5]|§| One could then extract the possible power-law exponent (as a function of 
\k — kpl) for fixed model parameters and compare to the expression in the second line 
of Eq. ([38D or Eq. (gO]) (the latter supplemented by Eq.(|T2D). 

In Refs. [19] and [2l] it was shown that G^^{x) (and also the stationary small 
momentum density correlation function) of the TL model can be computed as an average 
with a nonthermal statistical operator of a generalized Gibbs ensemble (GGE) with the 
eigenmode occupancies as the underlying set of integrals of motion. In the context of 
quenches the concept of GGEs was introduced to describe the (possible) stationary-state 
value of time evolved observables in systems with many conserved quantities [33l [Ml EH] . 
We verified that a similar GGE can be used if (i) the full momentum dependence of the 
potential is kept and (ii) for quenches between two interactions of arbitrary strength. 

3.2. A factor of two 

In systems which after a quench are expected to evolve into a thermal stationary state 
(described by a canonical ensemble), on intermediate time scales the appearance of pre- 
thermalized quasi stationary state was observed when computing the time evolution 
out of the noninteracting ground state in the weak coupling limit [261 122] • The quasi 
stationary state is characterized by observables oscillating for some time interval around 
a constant value which is different from the stationary one, the latter being reached for 
much larger times. Averaging observables like the fermionic momentum distribution 
function (at fixed k) over times in which the system is stuck in the pre-thermalized 
state gives values which agree to the corresponding ground state expectation values of 
the interacting system described by Hf up to a characteristic factor of two [261 EI]- 
It was later argued [28] that the pre-thermalized states correspond to the nonthermal 
steady state of systems with a sufficiently large number of integrals of motion. We 
thus expect to find these factors of two for the TL model, which falls into this class of 
systems, when considering the weak coupling limit. To keep this section compact we 
here exclusively consider the case of the quench out of the noninteracting ground state. 
Following Refs. [261 122] for the fermionic momentum distribution function we expect 

2[nUk)-n%{k)]=n^\k)-n%{k) (41) 

to hold up to order v'^. Here 

n% = e{kp - k) (42) 

§ One crucial difference of a generic lattice model to the TL model is the nonlinearity of the single- 
particle dispersion of the fermions. Similar to the momentum dependence of the two-particle potential 
this might have an effect on the transient dynamics. For a review on the effect of the nonlinear dispersion 
on equilibrium LL physics, see Ref. [32] ■ 
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is the noninteracting ground-state momentum distribution function (of right movers). 
Using the expansion (see Eq. ([9])) 

s\k) = ^v\k) + 0{v%k)) (43) 

in Eqs. ( IMl) and (l33l) we obtain 

2G-W-G-W = ^^ + 0(«r (44) 

The equihbrium ground-state Green function G'^'^(x) Eq. (IMl) is computed with the 
interaction after the quench that is with sf{k) — > s'^{k). Fourier transformation then 
gives Eq. ( |4T1) up to third order. 

We take the opportunity and also compare the anomalous dimensions of the steady 
state and the ground state characterizing the momentum distribution function close to 
kp. The former is given in Eq. (l38l) . the latter in Eq. ( l36l) . Expanding in the interaction 
using Eq. P3l) gives 

27gs -lst = {v^) . (45) 

Interestingly, we thus also obtain the relative factor of two in the power-law exponents 
of the momentum distribution functions at sufficiently weak coupling. 

Our comparison of the steady-state and ground-state momentum distribution 
function of the TL model provides additional evidence that the notion of a pre- 
thermalized state in a general model [261 [27] and its relation to the steady state in 
a model with many integrals of motion [28] is indeed meaningful. 

4. Time evolution after a quench 

We now evaluate the momentum integral of Eq. fl22l) for the Green function using the 
different two-particle potentials introduced at the end of Sec. [2] as well as with the ad 
hoc replacement Eq. ( !27|) (and the corresponding Eq. (l26|l ). In addition, we Fourier 
transform the Green function (see Eq. (fTTl) ). This way we obtain explicit results for the 
time evolution of the fermionic momentum distribution function after an interaction 
quench out of the noninteracting ground state. The results obtained for the different 
potentials are compared. 

4.I. Analytical insights 

4-1.1. Jump at kp We start out with an analytical result which can be obtained 
independent of the form of the two-particle potential (and even for the ad hoc 
procedure). Using Eq. (!22|) it is straightforward to show that the jump Z(t) of n{k,t) 
at kp defined as 



is given by 



Z{t) = lim n{k,t) - lim n{k,t) (46) 

t\i / nip n, XjfiJjT 

Z(t) = exp |- Fdk ^^ ^^^^ ^^^ (1 - cos [2w(A;)t])| . (47) 



Quench dynamics of the Tonionaga-Luttinger model 13 

For large times the remaining momentum integral can be performed using asymptotic 
analysis [31] leading to the result 

Z{t) ~ r^=' (48) 

first obtained for the ad hoc procedure in Ref . [19] . We here introduced the dimensionless 
time 

i = v^q,t. (49) 

4.1.2. The ad hoc procedure With the assumption Eq. (126|) and the replacement 
Eq. (1271) the momentum integral in Eq. (122]) can be performed analytically leading 

to 

Comparing Eq. fl32l) and the first line of Eq. fl38l) we find that g"^ = 7st and have 
accordingly replaced the former by the latter. Using the above expression the momentum 
distribution function can be written as 

An^<^^°^(A;,t) = n'^'^°^{k,t) - lim n'''^^°'{k,t) 

t—>oo 

iir r ,^ g'^'"'"^" [ [(^ + 2^Ft)^ + /g]""^' [{x - 2vptf + /g]""^' ^ 



2^ J~oo [X + 20] [X2 + /2]^-/2 I ^^2^^^^2 ^ ^2] ^-/2 

2'0FAA;t|/ 



sm 
dy 



IT (2{;Ft)^- io y [y2 + (2{;Ft)-2] 



7st/2 



^^^ / [(1 + y)^ + {2v,i)-T'^' [(1 - y)^ + (2i)pt)-^]"-/^ 1 
I [l + (2^pt>2]-=^/2 / • 

For 7st < 2 the integrand excluding the sine factor has a cusp at y = 1 — for t — > 00 the 
first derivative ioi y /^ 1 goes to minus infinity while for y \ 1 it goes to plus infinity. 
For fixed k ^ kp (for the behavior at k = kp, see above) this cusp leads to a nonanalytic 
power-law decay of An^'^^"'^ {k , t) given by 

An^'^^°=(A;,t) ~t-i-i^=' sin[2{)FAHJ , (52) 

which follows from restricting the integral to a small region around y = 1. The 
integral also has regular parts starting with a term ~ t~^. Therefore Eq. ( I52l) gives the 
asymptotic behavior as long as |7st < 1- Although it is possible to analyze An^'^^°'^{k, t) 
for stronger interactions we here refrain from doing so as we are mainly interested 
in the limit of weak to intermediate interactions. Loosely speaking for large t the 
smooth contributions away from the cusp average out due to the oscillatory term while 
close to the cusp the integrand changes quickly providing the leading contribution. 
This argument implies that the (dimensionless) time scale tp at which the power-law 
decay sets in increases for decreasing \k — kp] because the frequency of the ?/-oscillation 
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Figure 1. The time dependence of | An(fc, t)\ for the ad hoc regularization with v = 0.5 
at fixed AA; = 0.1 (red circles) and Ak = 0.01 (bhie crosses) on a double logarithmic 
scale. The solid lines show power-law fits to the envelope. Inset: Exponents extracted 
for diS'erent relative momenta Ak (circles: 0.8; squares: 0.4; triangles: 0.1) as a 
function of 7st, that is the interaction strength. The dashed line is the analytical 
result 1 + 37st/2. 



decreases: tp ~ AA; . The closer the momentum is to the nonanalyticity at kp the 
larger the time scale on which the asymptotic behavior sets in. In summary, within 
the ad hoc procedure the momentum distribution function at fixed k and for t ^ tp 
approaches its stationary value in an oscillatory fashion with (dimensionless) frequency 
2'OfA/c and an amplitude decaying as a power-law in t with exponent 1 + |7st (as long 
as the interaction is not too strong, that is for |7st < !)• 

These analytical insights can be confirmed numerically by performing the integral 
in Eq. (15T|) . Figure [1] shows | An'^'^'^°^(/i;, t)| for fixed AA; = 0.1 (red circles) as well as 
AA; = 0.01 (blue crosses) and v = 0.5 on a double-logarithmic scale. The solid lines 
are power-law fits (for times t G [400,4000]) to the envelope. In the inset exponents 
extracted along this line for different v (and thus 7st; see Eq. fl38|) ) and A A; are presented. 
They nicely fall onto the analytical prediction 1 + |7st shown as the dashed line. 
Consistent with the above analytical result the asymptotic behavior is reached faster 
the larger AA; (compare the data for AA; = 0.1 and 0.01 in the main panel of Fig. [1]). 



4.1.3. The box potential Analytical progress is also possible in the case of the box 
potential Eq. f l28|) . Then the argument of the exponential function in Eq. f l22|) simplifies 
to 

4s\ky{k) 



dk 



k 



(cos [kx] - 1) (1 - cos [2uj{k)t\) 



7st 







dk ^-^ (1 - cos [2vFkt\ ) 



(53) 
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Figure 2. The time dependence of |Ari(fc, t)\ for the box potential with u = 2 at fixed 
Afc = 0.4 (red circles). The solid line shows a power law ~ £~^ which fits the data for 
sufficiently large i ^ 200. Inset: Fourier spectrum of the data from the main plot. 



For asymptotically large t the integral can be evaluated using integration by parts 

1 7st 



/o k 2vFt 

From this it straightforwardly follows that 

sin (2vFt) 



(cos [qcx] - 1) sin {2vFt) + 0{t~ 



(54) 



An'^°"(A;,t) 



t 



(55) 



The two characteristic differences to the asymptotic dynamics of the ad hoc procedure 
Eq. ( 152|) namely the (i) independence of the decay exponent from the interaction strength 
and the (ii) independence of the oscillation frequency from the momentum k — kp aX 
which the distribution function is evaluated can be traced back to the nonanalyticity of 
the box potential; the long-time dynamics is dominated by the position of the jump in 
the potential, that is the upper boundary of the momentum integral in /. 

Figure |2] shows |Ar2'^°^(A;, t)| for fixed AA; = 0.4 and v = 2 obtained by numerically 
performing the momentum integral in Eq. (!22|) and the Fourier integral (with respect 
to position) Eq. ( ITTll . In the data a second smaller frequency than 2vf is observable. A 
Fourier analysis (with respect to time) shows that it is equal to the frequency found for 
^^adhoc^^^ t) and given by 2'OfA/c; see the inset of Fig. |2j In fact, the numerical data are 
consistent with the long-time dynamics for the box potential being dominated by the 
interplay of the term Eq. ( !55|) and another one of the form Eq. (!52|) . Considering a fixed 
time interval the t~^ decay is only clearly observable if the interaction and thus 7st is not 
too small and the term Eq. fl52|) can be neglected compared to Eq. fl55|) . Consistently, 
for such interactions the low-energy peak in the frequency spectrum carries a much lower 
weight than the high-energy one; see the inset of Fig. O 

It is possible to provide analytical evidence for the appearance of the second term 
Eq. ( I52l) at large t. To this end one uses that the argument / of the exponential function 
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in Eq. (122 p can be further evaluated 

/ = 7st [Ci{\q,x\) - 7 - ln(|g,a;|)] + %t [Ci(27)Ft) - 7 " H'^vpi)] 

- ^ [Cii\q,x + 2vFi\) - 7 - \n{\q,x + 2vFi\)] 

- ^ [Ci{\q,x - 2vFi\) - 7 - ln(|gcx - 2f;Ft|)] (56) 

with the integral cosine function Ci and the Euler constant 7. With this |Ar2^°^(A;, t)| 
can be brought into a form similar to Eq. ( ISTl) 



^ ^ vr (2{;Ft)7- Jo 2/^+^^ ^ ^ ' ' 



X 



exp {7st [Ci(2{;Ft) - Cii2vFi[y + l])/2 - Ci(2{;Ft|y - l|)/2] } - 1 • (57) 



Leaving out the sine factor and considering large t the integrand has similar to the 
one of Eq. (1511) a cusp at y = 1. This again leads to a nonanalytic term of the form 
Eq. (|52|) . In fact, the integrands (as a function of y) of the box potential and the ad hoc 
procedure coincide close to y = 1 up to the crucial difference that for the box potential 
the cusp-like behavior is modulated by an oscillation with frequency 2'OFt (associated to 
the oscillatory behavior of Ci). Taking everything together the appearance of the two 
terms Eqs. ( I52ll and ( l55l) is thus plausible from the analytics. 

To observe the power-law decay of \An^"^{k,t)\ for fixed k — kp in a given time 
interval (at large times) one again has to stay away from the nonanalyticity (the jump 



for finite t) at k = kp: the smaller 
power-law behavior sets in. 



Ak 



the longer it takes before the asymptotic 



4.1.4- Stationary points in u}{k) For the regular potentials Eqs. ( I29l) -( 13T]) the dispersion 
u}{k) Eq. (j9]) is a nonlinear function of the momentum. For each of the potentials a 
critical interaction strength Vc exists beyond which uj{k) has two stationary points. For 
the Gaussian potential one e.g. finds Vc = e^ ~ 7.39. For v > Vc the large time dynamics 
of Eq. (I22I) is dominated by these stationary points. Using the stationary phase method 
[3T] it is straightforward to show that on asymptotic time scales 

A fu +\ sin[2w(/i;i)t + 0i] sin [2w(/c2)t + ^2] .^^. 

An{k,t) r^ ci j= hC2 7= , (58) 

Vt Vt 

with amplitudes ci/2, phases 01/2, the two stationary points ki/2, and the dimensionless 
dispersion uj{k) = Co'(/c)/(t'FQ'c)llll Here we are primarily interested in the behavior 
at small to intermediate interactions (see also the above subsection on the ad hoc 
procedure) and thus focus on i) < Wc from now on. 

4.2. Numerical results 

For the Gaussian Eq. (129|) . the exponential Eq. ( 130|) . and the quartic potential Eq. (!3T|) 
we did not succeed in obtaining analytical results for An{k,t). The following analysis 

II For an example in which this type of behavior dominates for all interaction strengths, see Ref. |21| . 



Quench dynamics of the Tomonaga-Luttinger model 



17 



c 1.10 

CD 

o 

x1.05 

CD 



o gaussian 
X exponential 
A quartic 




Figure 3. Exponents of the power-law decay of \An{k,t)\ for different potentials 
(symbols) as a function of 7st, that is the interaction strength (see Eq. ([55)) ). The solid 
line shows 1 + 7st, the dashed one the exponent 1 + 37st/2 obtained within the ad hoc 
procedure, (a) Fixed momentum Afc = 0.8 and (b) Afc = 0.4. 



thus solely relies on the numerical evaluation of the momentum integral in Eq. ( |22l) and 
the successive Fourier integral Eq. f lT7|) . As the integrands are oscillatory functions one 
has to use routines which are adopted to this situation. Furthermore, for large times 
\An{k, t) I becomes very small (of order 10~^ and smaller) which requires a very accurate 
evaluation of the integrals. This limits the numerics and reliable results can be obtained 
up to times of the order of t = 5000. 

From the above analytical consideration we expect that (i) \An{k,t)\ decays as a 
power law in t and that (ii) this should be observable on moderately large times if Ak 
does not become too small. Our numerical results for |An(A;,t)| (not shown; the general 
form is similar to the data of Fig. [1]) are consistent with this expectation. In Fig. [3] we 
show how exponents extracted by a power-law fit of \An{k,t)\ for times t G [100,2400] 
depend on 7st, that is v (symbols). In Fig. |3]^a) the momentum is fixed at Ak = 0.8 
and in Fig. [3](b) at Ak = 0.4. The dashed line shows the result 1 + 37st/2 obtained for 
the exponent within the ad hoc procedure. Obviously the data for the three potentials 
coincide and differ from 1 + 37st/2 as well as the exponent 1 found for the box potential. 
In Fig. [3]^a) they instead nicely fall onto the line 1 + 7st. For smaller Ak = 0.4 the data 
slightly scatter around this line but are still consistent with it. The largest deviations 
are observed for the quartic potential (see below). Already at this stage of the analysis 
we can conclude that both the ad hoc procedure as well as the box potential fail in 
producing the exponent obtained for the more generic (regular) potentials. 

In Fig.m^a) we collected the data for the decay exponent of \An{k, t) | as a function of 
7st for the Gaussian potential from Figs. [3](a) and (b) and added another set obtained for 



A A; = 0.1 (open symbols). For small 



Ak 



the extracted exponents fall between the lines 



l+37st/2 and l+7st and one might be tempted to conclude that the exponent depends on 
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Figure 4. Exponents of the power-law decay of |An(fc,i)|. (a) The open symbols 
show the exponent as a function of 7st for the Gaussian potential and the momenta as 
given in the legend. The filled light blue triangles are for the quartic potential with 
Afc = 0.1. The solid line shows l + 7st the dashed one the exponent l + 37st/2. (b) The 
filled symbols show the dependence of the exponent on the upper boundary of the time 
interval over which the power law was fitted. The lower boundary is i = 200. The data 
were obtained for the Gaussian potential with v = 0.5, corresponding to 7st ~ 0.0417 
and Afc = 0.15, 0.3 as well as 0.4. 



k — kp. To further investigate this we studied the dependence of the extracted exponent 
on the time interval over which the power law was fitted. We increased the upper 
boundary t^ at fixed lower boundary ti. For the three regular potentials and all the Ak 
we studied we found that by increasing t^ the exponent tends towards l+7st- In Fig- 11(b) 
we show the t^ dependence of the extracted exponent for the Gaussian potential with 
ij = 0.5, corresponding to '-fst ~ 0.0417, and the relativ momenta Ak = 0.15, AA; = 0.3, 
and Ak = 0.4. For increasing fitting range with lower boundary ti = 200 the exponent 
has a clear tendency towards l+7st and appears to saturate for times we can reach within 
our numerics. We conclude that our results for the regular potentials are consistent with 
the assumption of an asymptotic power-law exponent which is independently oi k — kp 
given by 1 + 7st (as long as \k — kp\ 7^ 0). 

Figure Hl^a) additionally contains data for the exponent of the quartic potential at 
AA; = 0.1 (filled light blue triangles) extracted for the time interval t G [100, 2400]. The 
exponents are even closer to the result of the ad hoc procedure 1 + 37st/2 than the 
corresponding ones for the Gaussian potential at the same AA; (open green triangles). 
This can be understood as follows. First one realizes that it is the momentum 
dependence of the dispersion a;(A;) and not the one of the prefactor s^(A;)c^(A;) in Eq. (l22l) 
which dominates the long-time dynamics of the Green function. In the ad hoc procedure 
the nonlinear dispersion u}{q) is replaced by the linear one ufIq'I (see Eq. ( 1271) ). For small 
momenta the quartic potential deviates from its zero momentum value only to fourth 
order in q/qc- This implies that for small \q\/qc, ^{o) deviates from the linear behavior 
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to fourth order and over the momentum range < \q\/qc ~ 1 the dispersion relation of 
the ad hoc procedure and the quartic potential coincide. This has to be contrasted to 
the deviation from a linear dispersion for the exponential and Gaussian potentials which 
is of first and second order in q/qc, respectively. In this sense the quartic potential is the 
one considered closest to the ad hoc procedure|5ll For small Ak this implies, that an 
apparent power-law decay with exponent 1 + 37gt/2 is observable on intermediate times. 
Based on an analysis of the type discussed in the last paragraph we conclude that on 
asymptotic time scales the exponent will cross over to 1 + 7st. 

Similar to the cases of the ad hoc procedure and the box potential also for the 
more regular potentials An{k,t) oscillates around zero (with a decaying amplitude). A 
Fourier analysis (not shown) clearly reveals that the (dimensionless) frequency of the 
oscillation at fixed k — kp is given by 2u{k — kp). This is consistent with the result for 
the ad hoc procedure for which 2u{k — k-p) = 2vp\Ak\ because of the linearization. 

A detailed account of the time dependence of n{k, t) after a quench from one 
repulsive interaction to another one (obtained by Fourier transforming the Green 
function Eq. (!23|) ) can be obtained along the same lines. This is left for future work. 

5. Summary 

We have studied the quench dynamics of the TL model focusing on the single- 
particle Green function and the fermionic momentum distribution function. Instead 
of using an ad hoc procedure to regularize momentum integrals in the ultraviolet 
we kept the momentum dependence of the two-particle interaction and considered 
different potentials. The steady-state momentum distribution function close to kp 
is independently of the detailed form of the potential characterized by a power-law 
nonanalyticity with an exponent 7st which depends via the LL parameter K (or K^ 
and Kf when quenching between two repulsive interactions) on the strength of the 
interaction. We raised the question if this type of universality extends beyond the TL 
model to all models falling into the LL class in equilibrium. Importantly, 7st differs from 
the exponent 7gs characterizing the ground-state momentum distribution function at the 
same interaction. The ground-state and steady-state momentum distribution functions 
differ at weak coupling by a characteristic factor of two, known from systems in which 
pre-thermalized states appear. 

Our analytical as well as numerical results for the asymptotic time dependence of 
n{k,t) at small to intermediate interaction and k — kp ^ are consistent with 

An{k,t) L^. (59) 

For the class of generic regular potentials Eqs. fl29|) -( l3T|) the frequency i){k) is given by 
2uj{k — kp) and the power-law exponent by ,^ = 1 -|- 7st. While the frequency depends 

^ The box potential implies a strictly linear dispersion up to gc- The asymptotic dynamics still shows 
a different exponent compared to the one of the ad hoc procedure as the discontinuity of the potential 
strongly affects the long-time dynamics; see above. 
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on the details of the potential (the nonlinearity of the dispersion) the exponent of this 
class of potentials turns out to be universal in the LL sense; it is solely determined 
by the potential at zero momentum transfer and thus the LL parameter K. Both, the 
nonanalytic box potential as well as the commonly employed ad hoc procedure fail in 
reproducing this universal behavior. For the box potential ^ = 1 independent of the 
interaction strength and the (dimensionless) frequency is for all A/c given by 2{)f. In the 
ad hoc procedure we obtain a larger decay exponent .^ = 1 + 37st/2 and the frequency of 
the linearized dispersion i){k) = 2vp\Ak\. At k = kp, n{k, t) has a Fermi liquid-like jump 
of height Z{t) ~ t~'^=*. This holds independent of the potential considered. The time 
scale tp on which the power-law decay Eq. ( l59l) is observable at fixed k — kp increases 
if k approaches kp. If t ^ tp is not fulfilled the asymptotic behavior Eq. fl59l) might be 
superimposed or even hidden by additional contributions. 

In analogy to the posed question of LL-like universahty in the steady state it 
would again be very interesting to investigate if some form of universality in the time 
dependence extends beyond the TL model studying the quench dynamics of other models 
falling into the LL universality class in equilibrium. Two obvious candidates would be 
the decay exponent 1 + 7st of \An{k,t)\ found for regular potentials in the TL model 
and the exponent 7st of the jump at kp. 

Studying the TL model with momentum dependent two-particle potentials we here 
performed the necessary first step towards a more detailed understanding of universality 
in the quench dynamics of LLs. 
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